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Abstract 

We study the asymptotic scaling of the entanglement of a block 
of spins for the ground state of the one-dimensional quantum Ising 
model with transverse field. When the field is sufficiently strong, the 
entanglement grows at most logarithmically in the number of spins. 
The proof utilises a transformation to a model of classical probability 
called the continuum random-cluster model, and is based on a property 
of the latter model termed ratio weak-mixing. Our proof applies equally 
to a large class of disordered interactions. 



1 The quantum Ising model 

The quantum Ising model in a transverse magnetic field is one of the most 
famous examples of exactly solvable one-dimensional quantum models. The 
solution was first given by Pfeuty in (26j, based on earlier works by Lieb, 
Schultz, and Mattis [18] and by McCoy [21J. The diagonalisation of the Hamil- 
tonian and the determination of the energy eigenstates is based on methods 
developed by Jordan and Wigner [TH| in the theory of second quantisation of 
fermion fields, and by Bogoliubov [7] in the theory of superconductivity. This 
model exhibits a second-order phase transition in the ground state when the 
temperature of the system is zero. The existence of the phase transition and 
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the computation of the spin-spin correlation functions were studied in |26j : 
rigorous results for the correlation functions in the presence of disorder are 
provided in [H [9] . 

Quantum systems, unlike classical systems, can support composite pure 
states for which it is impossible to assign a definite state to two or more sub- 
systems. States with this property are known as entangled states and have 
attracted a great deal of interest recently due to their resource-like proper- 
ties. The investigation of the entanglement properties of strongly interacting 
quantum spin systems, with a view toward quantum phase transitions, was 
initiated by Osterloh et al. [25] and by Osborne and Nielsen [23] (see, for 
example, [4] and the references therein for further studies). It is now under- 
stood that the strength of quantum entanglement is related to the number 
of parameters required to describe a quantum state classically. Thus, for ID 
systems, the scaling of the geometric entropy — the degree of entanglement 
of a distinguished subsystem with respect to the rest — has emerged as the 
crucial parameter which quantifies whether the state is hard or easy to sim- 
ulate [2H]. It has been conjectured that the entropy of entanglement obeys 
an area law, scaling as the boundary area in the subcritical phase, with a 
possible logarithmic correction for the critical phase. There is a paucity of 
rigorous results concerning the scaling of the entanglement of a block for the 
quantum Ising model; the above results are typically obtained by numerical 
calculations, or conformal field theory methods [3|. There are some rigorous 
derivations of the scaling of the entropy function for certain ID spin models 
(specialised essentially to the XY model), see [I] for further references. 

In this paper, we utilise a new method for studying the entanglement 
properties of the quantum Ising model. This is based on a representation 
formulated by Aizenman, Klein, and Newman [1] of the model in terms of a 
continuum random-cluster model on a certain space-time graph. (See also the 
earlier paper [9].) Using a technique termed ratio weak-mixing, developed by 
Alexander [2j [3] for random-cluster and Potts models on discrete lattices, we 
prove a bound on the entanglement entropy in the subcritical regime, when 
the magnetic field intensity is strong compared to the spin coupling. 

The quantum Ising model is defined as follows. Let L > 0. For m > 0, let 
A m = {— rn, —m + 1, . . . , m + L} be a subset of the one- dimensional lattice Z, 
and attach to each vertex x G A m a quantum spin-| with local Hilbert space 
C 2 . The Hilbert space 7i for the system is H = <S)™=- m < ^ 2 - A convenient 

basis for each spin is provided by the two eigenstates |+) = ( J , | — ) = 
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of the Pauli operator 

1 
-1 



a (3) 



at the site x, corresponding to the eigenvalues ±1. The other two Pauli 
operators with respect to this basis are represented by the matrices 

° I V a?) = ( ° 'J 



1 J ' x \ i 

A complete basis for Ti is given by the tensor products (over x) of the eigen- 
states of (j x ■ In the following, \(p) denotes a vector and (<f>\ its adjoint. As a 
notational convenience in this paper, we shall represent sub-intervals of Z as 
real intervals, writing for example A m = [—m, m + L] . 

The spins in A m interact via the quantum Ising Hamiltonian 

(x,y) x 

generating the operator e~^ Hm where f3 denotes inverse temperature. Here, 
X Xy y > and 8 X > are the spin-coupling and external-field intensities, respec- 
tively, and J2( xy ) denotes a sum over all (distinct) unordered pairs of spins. 
We concentrate here on the case of interactions between neighbouring spins: 
\ XtV = for \x — y\ > 2. While we shall phrase our results for the translation- 
invariant case X XtX +i = A and S x = 5, our approach can be extended to random 
couplings satisfying the condition 

F(X XtV < A) = > 6) = 1, (1.3) 

with 9 = X/5 a sufficiently small constant (see Section [H]). The ensuing Hamil- 
tonian has a unique pure ground state \ip m ) defined at T = (/? — > oo) as the 
eigenvector corresponding to the lowest eigenvalue of H m . In the translation- 
invariant case the ground state \ip m ) depends only on the ratio 9. 

For definiteness, we shall work here with a free boundary condition on A m , 
but we note that the same methods are valid with a periodic (or wired) bound- 
ary condition, in which A m is embedded on a circle. One difference worthy 
of note is that the correlation functions of the critical model are expected to 
depend on the choice of boundary conditions, see [26] . 

We write p m (0) = e~ pHm / tr(e _/3Hm ), and 

p m = lim p m (0) = \ip m ){ip m \ 

for the density operator corresponding to the ground state of the system. The 
existence of the limit follows by random-cluster methods, see [I| , and we return 
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to this in Section [31 The ground-state entanglement of \ip m ) is quantified by 
partitioning the spin chain A m into two disjoint sets [0, L] and A m \ [0, L] and 
by considering the entropy of the reduced density operator 

p L m = tr Am \[o,L](|<M(^m|)- (1-4) 

One may similarly define, for finite /3, the reduced operator In both 

cases, the trace is performed over the Hilbert space ((£>~=_ m C 2 )<E>((S)^J!+i C 2 ) 
of the spins belonging to A m \ [0,L]. Note that is a positive semi-definite 
operator on the Hilbert space Tii of dimension d = 2 L+1 of spins indexed 
by the interval [0, L]. By the spectral theorem for normal matrices [5J, this 
operator may be diagonalised and has real, non-negative eigenvalues, which 
we denote Aj(p£j. The arrow indicates that the eigenvalues are arranged in 
decreasing order. 

Definition 1.5. The entanglement of the interval [0, L\ relative to its com- 
plement A m \ [0, L] is given by 

S{p L m ) = -ti{ P L m \og 2P L m ). (1.6) 

This quantity may be expressed thus in terms of the eigenvalues of p^: 

2 L+1 

S(P L J = -£ Aj(p^) log 2 Aj(^), (1.7) 
j=i 

where log 2 is interpreted as 0. 

In Section El we prove our main theorem: the order of the entanglement 
scaling is at most log 2 L for the ground state in the sub critical regime. This 
result follows as a corollary of the main estimate, given by Theorem 16.51 in 
Section [61 In Sections [3H11 we describe the mapping of the density operator of 
the quantum Ising model to a stochastic integral in terms of a Poisson measure 
(as in P). The mapping begins by considering states with (3 < oo and deriving 
the ground state in the limit (3 — ► oo. This allows us to express the matrix 
elements of the ground state in terms of a continuous percolation model on a 
two-dimensional space-time graph, with one continuous axis describing time. 
In this setting, the elements of the reduced state are related to a random- 
cluster model on the same graph, but with the addition of a 'slit' along the 
interval [0, L] at time 0. The continuum random-cluster model is presented 
in detail in Section [5) Section [6] contains the main result, which allows us 
to establish the scaling of the entanglement entropy, while in Section [7] we 
explain the technique of ratio-weak mixing on which the proof is based. The 
extension of our results to disordered systems is discussed in Section [HI 
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2 Entropy of the reduced state 

In this section, we study the behaviour of the entropy of the reduced state p^ 
in the subcritical regime (with 9 = X/S small). In order to derive an adequate 
upper bound on the entropy, we shall analyze the influence on the spectrum of 
the reduced density operator produced by imposing a change in the boundary 
conditions of the spin chain. Specifically, we consider the distance between the 
largest eigenvalues of two states defined on [0, L] with respect to two different 
lattices, A m , A n , with m <n. The entropy will be estimated by studying the 
operator norm 

\\p L m -p L n \\ = sup \(i'\pi- P^)\, (2.1) 
IHI=i 

where the supremum is taken over all vectors G Hl with unit L 2 -norm 
belonging to the Hilbert space Hl of spins in [0,L]. We shall see in Section 
H] that \\pm — p%\\ may be expressed in terms of a certain random-cluster 
representation of the quantum Ising model. In Sections [6] and [7] we shall use 
a coupling of random-cluster measures and the method of 'ratio weak-mixing' 
to prove the following. 

Theorem 2.2. Let A, 5 G (0, oo) and write 9 — X/8. There exist constants 
a,C G (0, oo) depending on 9 only, and a constant 7 = 7(6') satisfying < 
7 < 00 if 9 < 1, such that, for all L > 1, 

\\p L m -p L n \\ < min{2,OL a e-^}, 2 < m < n. (2.3) 

Furthermore, we may find such 7 satisfying 7 — > 00 as 9 j 0. 

Proof. That p%\\ < 2 is a consequence of the fact that the are density 

operators. An upper bound of the form C'L a e~ im holds by Theorem 16.51 and 
the preceding discussion whenever m > M for suitable M = M(9). Inequality 
(TO) follows on replacing C by C = e 7M max{C", 2}. □ 

We shall apply (I2.3P iteratively in order to obtain an upper bound for the 
decay of the vector of eigenvalues {Aj(p^) : j = 1,2,...}, valid for all large 
m. The proof makes use of the following decomposition property, valid for 
any pure state of a bipartite system, see [23] . 

Theorem 2.4 (Schmidt decomposition). Let \ip m ) be the pure ground state 
of the composite system [0, L] U (A m \ [0, L]). There exist orthonormal bases 
{\ u j){o,L]i \ v k) A m \{o,L]} f or the states of [0, L], A m \ [0, L] respectively, such that 

IVV> = yfiiOk) K)[0,L]NA m \[0,L], (2.5) 

where s, the Schmidt rank, is given by s = min {2 L+1 ,2 2m }. 
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Proof. We begin by writing \if) m ) in terms of an orthonormal basis \a) [o,l] |/?) A m \[o,L] 
where |a)[o,L] (respectively, \/3) A m \[o,L]) is an orthonormal basis for the spins 
in [0,L] (respectively, A m \ [0, L]): 

a=l (3=1 

where 

2-Zj + I 2^ m 
Q=l 0=1 

The coefficients constitute a 2 L+1 x 2 2m matrix and, as such, we can 
apply the singular- value decomposition [6] to write 

s 

3=1 

where s = min{2 L+1 , 2 2m }, U a j is a 2 i+1 x s-sized isometry, dj > for all 
j = 1,2,..., 2 L+1 , and Vjp is an s x 2 2m -sized isometry. Defining 

s s 

\ U j)[0,L] = s]Uqj\a)[0,L], |^)a ot \[0,L] = ^2Vjp\P)A m \[0,L], 
a=l f3=l 

we see, because U and V are isometries, that {|wj)[o,l]} an d {\ v j) A m \[o,z]} are 
orthonormal sets of vectors for the spins in [0, L] and A m \ [0, L], respectively. 

A simple computation shows that the reduced density operator for the 
spins in [0, L] is given by 

s 

P™ = E4KW]<^I (2-6) 

3=1 

and so we identify dj = sj Aj(p^), after re-ordering the index j if necessary. 
Note that the rank of is less than or equal to the Schmidt rank of \ip m ). □ 

We compute the entanglement of [0, L] with respect to the rest of the 
system as in (11.71) . 

s 

S{p L m ) = -£ \\{p L m ) log 2 \\{p L m ). (2.7) 

3=1 

Here is our main theorem. With the exception of the natural logarithm func- 
tion In, all logarithms in the remainder of this section are taken to base 2. 
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Theorem 2.8. Consider the quantum Ising model ( 11.21) on n = 2m + L + 1 

spins, with parameters X, 5, and letj, a, C be as in Theorem \2.^A If'y > 41n2 ; 
there exist constants C\ and c 2 depending on 7 only such that 



In summary, the entanglement entropy S(p^) is at most logarithmic in L 
if the field strength 5 is sufficiently large. The bound 4 In 2 is sufficient but 
not necessary, and may be improved with more care in the proof. We do not 
know how to replace this condition by 7 > 0. 

We believe that the upper bound (12.91) is, in many cases, not tight. For 
the translation-invariant subcritical case 8 = X/5 < 2 it is expected, on phys- 
ical grounds, that the upper bound can be improved to a constant. (See [1] 
and the references therein for an extensive review of the physical arguments 
for entanglement scaling in non-critical and critical quantum spin models.) 
Renormalisation group arguments and conformal field theory methods sug- 
gest that, at a critical point, the upper bound should scale with logL. For 
9 > 2 the system enters the supercritical regime where the system has two 
degenerate ground states, and the ground state is no longer a pure state. 
Nonetheless, it is expected that the entropy of a block is again bounded by 
a constant. For higher dimensions d > 2 our argument breaks down because 
the number of non-zero Schmidt coefficients for a distinguished region grows 
too quickly for our perturbation argument. 

The proof of Theorem 12.81 follows an iterative inductive procedure, where 
at each step the distance k from the boundary of [0, L\ is increased and the 
spectrum of the relative density operator p\ is estimated. We illustrate the 
procedure by the following simple case: consider the ground state | V^o) for the 
Ising model defined on only L + 1 spins. In this case the reduced density 
operator pg for [0, L] is exactly Pq = \4>o)(il>o\, i.e., a pure state, with entropy 
S(pq ) = 0. When m — 1, the reduced density operator p\ for the region [0, L] 
is mixed, but it has at most 2 2 non-zero eigenvalues. This follows from the 
Schmidt decomposition applied to the ground state across the bipartition 
[0, L] U (Ai \ [0, L]). Thus, the entropy of the block [0, L] is bounded above by 
S(p\) < 2. Consider now the reduced density operator p%. By the Schmidt 
decomposition, the operator p£ has at most 2 2k non-zero eigenvalues. Assume 
that 2k < L + l, and consider the addition of a single spin at either boundary. 
The new reduced density operator p£ +1 has at most four times as many non- 
zero eigenvalues as p£. However, by (12.31) . 



S(pt) < cilog 2 L + c 2 , 



m > 0. 



(2.9) 



so that the eigenvalues of p\ 




(2.10) 
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Proof of Theorem \2.8i . Let K = [7" 1 ln(CL a )] , with C, a, 7 as in Theorem 
12.21 We shall assume that m, K > 2, 7 > 4 In 2. There are two cases, depend- 
ing on whether m < K or m > K. Assume first that 2 < m < K. The rank 
of equals the Schmidt rank 2 2m of \ip m ). Therefore, 

{ 2 2m N 

where the supremum is over all non-negative sequences p = (pj : 1 < j < 2 2m ) 
with sum 1. Hence, 

S(ffa) < log 2 2m = 2m < 2K, m < K. (2.11) 

Assume next that m > K. We shall apply the following theorem, see [6]. 

Theorem 2.12 (Weyl perturbation theorem). For Hermitian operators A and 
B on a Hilbert space of dimension n, 

max|Aj(A)-Aj(£)| < \\A-B\\. (2.13) 

Let e(r) = CL a e~ 1 ^ K+r \ and note by the definition of K that 

e(r) < e _7r , r > 0. (2.14) 
Setting A = p L K , B = p% +l in Theorem |2H2l we deduce by (EOl that 

max|Aj(p|)-Aj(pi +1 )| < e(0). (2.15) 

Therefore, 

|Aj(^ +1 )|<Aj(pi) + e(0), j = l,2,...,2 2 *, 

|AjG4 +1 )| < 6(0), j = 2 2K + 1, 2 2 ^ + 2, . . . , 2 2 ^ +1 \ (2.16) 

We shall now iterate this process in order to obtain a bound on the eigen- 
values of Px +r , for r > 1. There are three cases: 

(i) j < 2 2K , in which case 

r-l 

Ai(p£ +r )<Aj(p£) + 5>(0; (2.17) 

1=0 
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(ii) 2 2K < 2^ K+S ^ <j< 2 2 ^ +s+1 ) < 2 2 ^ K+r \ in which case 

r-l 

X l M <+r )< (2.18) 

(iii) 2 2 ^ +r - ) < j, in which case 

Aj(pi- +r .) = 0. (2.19) 
Let s = [\ log j'J - if, so that, by QgHj) , 

oo 

A)(p^)<Aj(pi) + ^e-^, j<2 2 *, 
z=o 

oo 

which is to say that 

Aj(p^)<Aj(p^) + c , j<2 2 *, 

Aj(p, L J < c e-^, 2 2 *<j, (2.20) 

where 

1 - e-T ~ 3 



14 , 
Co = , < 77- (2.21) 



By (E2D, 

Aj(p£) 2 2 *<j, (2.22) 

where £ = 7/(2 In 2) > 2 and c' Q = c' (L) = c e^ K+1) . 
By dLTD, 

S(p£j = S x + S 2 , (2.23) 

where 



1/ 2 L +! 



where v (> 2 2( - K+2 >) is an integer to be chosen later. We shall bound Si and 
S2 separately. Since the Aj(p^), 1 < j < v, are non- negative with sum Q 
satisfying Q < 1, 

Si < log 1/. (2.24) 



10 



G. R. Grimmett, T. J. Osborne, P. F. Scudo 



We shall use the tail estimate ( 12. 22ft to bound S 2 , making use of the fact 
that the function f(x) = — 2 log 2 satisfies: /(0) = 0, and f(x) < f(y) when- 
ever < x < y < e -1 . 

By (l2Tm and the definition of f , 



and so, recalling that v > 2 2( - K+2 ^ and £ > 2, 

s^-E4iog(|)<-f:^° g (f 

j=i/+i v 7 j=H-i v 

00 1 ^. , 00 , 

^ r ' i 'i 1 £ c o 

<-[cologc ] L^Lj 

< IcnlogCnl / — rdi + : — ^ / —\nxdx 
J u xt ln2 y v a* 



< 



We now set v = 

[ e 7(^+i)] to obtain 

S 2 < c x K + c 2 , (2.25) 

for suitable constants c±, c 2 depending on 7 only. By (I2.23H ( T2.25[) . 

5(p£)<ciif + 4 m>K, (2.26) 

which may be combined with (12.111) to obtain (12.91) with adjusted constants. 

□ 



3 Percolation representation of the ground state 

Aizenman, Klein, and Newman pQ derived a random-cluster representation 
for the thermal state of the quantum Ising Hamiltonian ( 11.21) . thereby relating 
spin-correlation properties to graph-connectivity properties. In this represen- 
tation, the thermal density operator, defined as 



Entanglement in the quantum Ising model 



11 



is described by a stochastic integral with respect to a Poisson measure. This 
Poisson measure is defined on the space-time graph A mj( g = A m x [0,/3], gen- 
erated by associating a continuous (imaginary) time variable t G [0, 0\ to each 
site x e A m . We refer to a line of the form {x} x [0, j3] as the time-line at the 
site x. 

For completeness, we reproduce here the derivation of the random-cluster 
representation of the ground state, and we derive the corresponding represen- 
tation for the reduced state on [0, L\. Note that the derivations are valid with 
the line A m replaced by any finite graph G. By ( II. 2p with v — \ ^2/ xy \ AI and 
I the identity operator, 

e -W m +») = e -W+v) } (3.2) 

where 

l7 = -*2>« ^=-i^A(^ 3 )-I), 

x {x,v) 
and the second summation is over all neighbouring pairs in A m . Although 
these two terms do not commute, we may use the so-called Lie-Trotter formula 
(see, for example, [27]) to factorize the exponential in (jl.2p into single-site 
and two-site contributions due to U and V, respectively. By the Lie-Trotter 
formula, 

e - {U+ V)At = e -UAt e -VAt + 0(A ^ 

We divide the interval [0, /3] into N parts each of length At = 1/N, and deduce 
that 

e -P(P+v) = lim f e -UAt e -VAnl3/ A \ (3.3) 
We then expand the exponential, neglecting all terms of order o(At), to obtain 

e -(3(H m +u) 

/ v /3/A* 

&, ( n [c 1 - 5At ) n + 5Atp -] n K 1 - AAt ) n + AAtp iJ ) ' ^ 3 - 4 ) 

V 31 {x,y) ) 



where P, 1 = a\ x) + I and P 3 , = ±(a®a® + I). 

Let B be the set of basis vectors \o~) for ?i of the form \o~) = §Z) X |±) x - 
There is a natural one-one correspondence between B and the space P = 
nC^-^i - 1) +!}• We shall sometimes speak of members of P as basis vectors, 
and of H as the Hilbert space generated by P. Similarly, the space Hl of 
spins indexed by the interval [0, L\ may be viewed as being generated by 

ft = itSu{-i.+i>- 
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The stochastic-integral representation may be obtained from (13.41) by in- 
serting the resolution of the identity 

J>)H = I (3.5) 

aeP 

between any two factors of the products. The product (13 .4p contains a col- 
lection of operators acting on sites x and on neighbouring pairs (x, y) . By la- 
belling the time-segments as Ati, At 2 , • • • , At at in [0, j3], and neglecting terms 
of order o(At), we may see that each given time-segment arising in (13.41) con- 
tains one of: the identity I; an operator of the form P^; an operator of the form 
P xy . Each such operator occurs in the time-segment with a certain weight. 

Let us consider the action of these operators on the states |<r) for each 
infinitesimal time interval Atj, i £ {1 ; 2, . . . , N}. The matrix elements of each 
of the single-site operators are given by 

(v'\a^ + I\a)=5 a ,, ax +5 a ,, Wx = l, (3.6) 

where a x is the value of the spin at x in the (product) basis vector |<t), and 
W x is the opposite spin to a x . When it occurs in some time-segment Ati, we 
place a mark in the interval {x} x Ati, an d we call this mark a death. Such a 
death has a corresponding weight 5 At + o(At). 

The matrix elements involving neighbouring pairs (x, y) yield 

l(ay y \a x Va^ + l\a x a y ) = 8^,8^,8^. (3.7) 

When this occurs in some time-segment Ati, we place a connection, called a 
bridge, between the intervals {x} x Ati and {y} x Ati. Such a bridge has a 
corresponding weight A At + o(At). 

In the limit At — > 0, the spin operators generate thus a Poisson process 
with intensity 5 of deaths in each time-line {x} x [0, j3], and a Poisson process 
with intensity A of bridges between each pair {x} x [0, (3], {y} x [0, j3\ of time- 
lines, for neighbouring x and y. This is an independent family of Poisson 
processes. We write D x for the set of deaths at the site x, and B x>y for the set 
of bridges between neighbouring sites x and y. The configuration space is the 
set flm.p containing all finite sets of deaths and bridges, and we may assume 
without loss of generality that no death is the endpoint of any bridge. 

For two point (x,s),(y,t) G A m ^, we write (x,s) <-> (y, t) if there ex- 
ists a path from the first to the second that traverses time-lines and bridges 
but crosses no death. A cluster is a maximal subset C of A m p such that 
(x, s) <-» (y, t) for all (x, s), (y, t) G C. Thus the connection relation <-» gener- 
ates a percolation process on A = A m /3 , and we write Pa 5 a,5 for the probability 



Entanglement in the quantum Ising model 



13 



measure corresponding to the weight function on the configuration space Vt m ^. 
That is, Pa.a,<5 is the measure governing a family of independent Poisson pro- 
cesses of deaths (with intensity 8) and of bridges (with intensity A). The 
ensuing percolation process has been studied in [5]. 

We shall later need to count the number of clusters of a configuration 
uj G Q m ,p subject to any of four possible boundary conditions, of which we 
specify two next (the other two appear in the next section). The meaning 
of periodic boundary condition is that any clusters containing two points of 
the form (x, 0) and (x,/3), for some x G A m , are deemed to be the same 
cluster, and they contribute only 1 to the total cluster count. The meaning 
of wired boundary condition is that any clusters containing two points of the 
form (x, 0) and (y,(3), for x,y G A m , are deemed to be the same cluster and 
contribute only 1 to the total count. We write k v {uj) (respectively, k w (uj)) 
for the number of clusters of uj subject to the periodic (respectively, wired) 
boundary condition. Note that k w (uj) — 1 is the number of clusters of uj (with 
free boundary conditions) that do not intersect [—m,m + L\ x {0, /?}. 

Equations (13.6| - fl3.7p are to be interpreted as saying the following. In 
calculating the operator e~ /3 ^ Hm+u \ one averages over contributions from re- 
alizations of the Poisson processes, on the basis that the quantum spins are 
constant on every cluster of the corresponding percolation process, and each 
such spin-function is equiprobable. 

More explicitly, 

e-^» + ") = I d¥ AXS (uj) It J] J] P^Kvit')) > (3-8) 

\ {x,t)dD ({x,y),t')eB J 

where T denotes the time-ordering of the terms in the products, and B (re- 
spectively, D) is the set of all bridges (respectively, deaths) of the configuration 
& £ ^m,/3- The P^it) and P^ y (t) are to be interpreted as the relevant operators 
encountered at the deaths and bridges of uj. 

Let uj G VL m fi. Let E(a>) = S m i (u;) be the space of all functions s : 
A m x [0,/3] — > { — 1,+1} that are constant on the clusters of uj, and let 
be the counting measure on E(o>). Let K(uj) be the time-ordered product of 
operators in (13.81) . We may evaluate the matrix elements of K(uj) by inserting 
the resolution of the identity between any two factors in the product, obtaining 
thus that 

(a'\K(uj)\a) = Hs(;0) = <t}1{s{;P) = o'}, a, a' G P, (3.9) 

where 1{^4}, and later 1^, denotes the indicator function of A. This is the 



14 



G. R. Grimmett, T. J. Osborne, P. F. Scudo 



number of spin-allocations to the clusters of lu with given spin-vectors at times 
and p. 

The matrix elements of the density operator p m {P) are therefore given by 
(JlpmtfM = 4~ I Hs(;0) = a}l{s(;(3) = a'}d^(s)dF A , Xi5 (Lu), (3.10) 

An J 

for a, a' G P, where 

Z m = Z m (J3) = tr(e-^ + ")) (3.11) 
is the partition function. Thus, 

(o-'\ Pm {(3)\o-) = ±- fdP AtX>s (uj) J2 l{s(;0) = v}l{o-(;(3) = a'} 

m seE(w) 

= 4- f d¥ AX5 (uj) 2 feW ^- 1 l£ (CTiCT , ) H, a,a' G P, (3.12) 

where the final term in the integrand is the indicator function of the event 
E(a, a') containing all uj G £L m $ such that: for all x, y G [— m, m + L\. 

(x, 0) (y, 0) whenever a x ^ a y , 
(x,P) ^ (y,p) whenever a' x ^ a' y , 
(x,0) (y,P) whenever a x ^ a' . 

See Figure [1] for an illustration of the space-time configurations contributing 
to the Poisson integral (13.121) for the matrix elements of p m (P). 

On setting a = a' in (13.121) and summing over a G P, we find that 

Z m = tr(e-^ H - + ^) = [ 2 kP ^ d¥ AX5 (uj). (3.13) 



A method was developed in [T] (as amplified in the next section) to rep- 
resent (a'\p m (P)\a) as a certain probability, and to prove that it converges as 
P — > oo. In particular, it was shown in pQ that the ground state p m = \ip m ) ("0ml 
satisfies 

p m = lim — e -ft Hm+u \ (3.14) 

/3^oo Z m 



4 Percolation representation of the reduced 
state 

The analysis of the last section may be repeated for the reduced density opera- 
tor Pm(P) by tracing ( 13. 8ft over a complete set of states of the spins indexed by 
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Figure 1: An example of a space-time configuration contributing to the Pois- 
son integral (13. 12[) . The cuts are shown as circles and the distinct connected 
clusters (each of which contributes a factor 2 to the term 2 fcW( ^) are indicated 
with different line-types. 



A m \ [0, L]. The corresponding boundary condition for the configuration uj G 
Q m .f3 turns out to be partially periodic, in that any two clusters of uj containing 
points of the form (x, 0) and (x, for some x G [—m, — 1] U [m + 1, m + L], are 
deemed to be the same cluster and contribute only 1 to the total cluster count. 
No such assumption is made for sites x G [0, L], and we refer to the boundary 
condition on [0, L] as free. Let k pp (uj) be the number of clusters of uj subject 
to the partially periodic boundary condition. We shall need a fourth way to 
count clusters also, as follows. The periodic / wired boundary condition is that 
derived from the partially periodic condition by the additional assumption of 
a wired condition on [0, L\. any two clusters of uj containing points of the form 
(x, 0) and (y, (3), for some x, y G [0, L], are deemed to be the same cluster and 
contribute only 1 to the total cluster count. We write k pw (uj) for the number 
of clusters with the periodic/wired boundary condition. Note that k pw (uj) — 1 
is the number of clusters of uj (with the partially periodic boundary condition) 
that do not intersect [0, L\ x {0, /?}. 

As in (I3.10p - (13.12p . the matrix elements of the reduced state Pm(P) are 
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spin configuration cr 




spin configuration a' 



Figure 2: An example of a space-time configuration contributing to the ma- 
trix elements for the reduced density operator p^{P). The box has partially 
periodic boundary conditions and is drawn in such a way that the slit Si is 
at the centre. The spin configurations on the top and the bottom of the cut, 
and the connected clusters for this new cut geometry are indicated. 



given by 

(<y' L \p L MM = ±- I dF A>x Au)2 kPW ^-H EM) (u), a L ,a' L G P L , (4.1) 

where E(<j L ,a' L ) is the event that: if x, y G [0, L] are such that a^^ x ^ cr' L y 
then (x, 0) <^ pp (y,{3)- Here, <-> pp denotes the connectivity relation subject to 
the partially periodic boundary condition. See Figure [2] for an illustration of 
the slit space-time, and of the connected clusters contributing to the matrix 
elements of 

We shall study the entropy of the reduced state via the operator norm of 
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( 12.11) . Let \ip) G T~Ll have unit L 2 -norm, so that 

1^) = E C ( a L)(TL 

(r L eP L 



for some function c : Pl — > C with ^ ul&Pl c{gl)c{gl) = 1. Then 



(4>\f&(0)\1>) = E Wi)^ °i) (4-2) 



where 



pfaM = -^K\pim<J L ), a L ,a' L G P L , (4.3) 
*m= £ KIP™(/5)kL), (4.4) 



and 



a mij g = Z m /N m . (4.5) 

We shall see in the next sections that ( 14. 2D — (14. 41) may be written in terms of 
a certain probability measure on fl m> p called the random- cluster measure. 



5 The continuum random-cluster model 

Perhaps the best way to express the percolation representations of the ground 
and reduced states is in terms of the so-called random-cluster model on Z x R. 
We summarise the definition and basic properties of this model in this section, 
using the language of probability theory. The remaining part of the paper is 
a self-contained account of the model, and includes the proof of Theorem 12.21 
see Theorem 16.51 Of special interest will be the property of so-called ratio 
weak-mixing, studied earlier for the lattice case in [21 E]. 

We shall consider the (two-dimensional) continuum random-cluster model 
on the 'space-time' subset ZxRof the plane. The underlying space is {(x, t) : 
x G Z, t G R}, and we refer to Z as the space-line and R as the time-line. 
Everything proved here has a counterpart, subject to minor changes, in the 
more general setting of 7L d x R with d > 2, but we shall restrict ourselves to 
the case d = 1. 

We shall construct a family of probabilistic models on Z x R. Let A, 5 G 
(0,oo). In the simplest such model, we construct 'deaths' and 'bridges' as 
follows. For each x G Z, we select a Poisson process D x of points in {x} x R 
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with intensity 5; the processes {D x : x G Z} are independent, and the points 
in the D x are termed 'deaths'. For each x G Z, we select a Poisson process 
£> x . of points in {x + |} xl with intensity A; the processes : x G Z} 
are independent of each other and of the D y . For each x G Z and each 
(x + |,t) G -B x , we draw a unit line-segment in R 2 with endpoints (a;, t) and 
(x + l,t), and we refer to this as a 'bridge' joining its two endpoints. For 
(x,s), (y,t) G Z x R, we write (x, s) «-> (?/,£) if there exists a path tt in R 2 
with endpoints (x,s), (y,t) such that: 7r comprises sub-intervals of Z x 1 
containing no deaths, together possibly with bridges. For A, A C Z x R, we 
write A <-> A if there exist a G A and 6 G A such that a <-> 6. 

For (x, s) 6 Z x 1, let be the set of all points (y, i) such that (x, s) <-» 
(y, i). The clusters C XiS have been studied in [5], where it was shown in 
particular that 

/, , s f = 1 if < 1, , , 

P A5 C <oo - 5.1 

' \< 1 if 9 > 1, 

where = (0,0) is the origin of Z x R, 6* = A/5, and |C| denotes the (one- 
dimensional) Lebesgue measure of the cluster C. The process thus constructed 
is a continuum percolation model in two dimensions. As noted in [5], it differs 
from the contact model on Z only in that two points may be joined in the 
direction of either increasing or decreasing time. See [T91 12U] for details of the 
contact model. 

Just as the percolation model on a lattice may be generalised to the so- 
called random-cluster model (see [H]), so may the continuum percolation 
model be extended to a continuum random-cluster model. We shall work here 
mostly on a bounded box rather than the whole space Zxl, Let a, b G 
Z, s, t G M satisfy a < b, s < t, and write A = [a, b] x [s, t] for the box 
{a, a + 1, . . . , b} x [s,t] of Z x R. Its boundary <9A is the set of all points 
(x, y) G A such that: either x G {a, 6}, or y G {s, £}, or both. As sample space 
we take the set comprising all finite subsets (of A) of deaths and bridges, 
and we assume that no death is the endpoint of any bridge. For uj G f^A; we 
write B{ui) and D{u) for the sets of bridges and deaths, respectively, of uj. We 
take as ex-field Tk that generated by the open sets in the associated Skorohod 
topology, see [5], [10] . 

In order to maintain the link to the quantum Ising model, we choose to 
impose a top/bottom periodic boundary condition on A; that is, for every 
x G [a, b], we identify the two points (x, s) and (x, t). The remaining boundary 
of A, denoted <9 h A, is the set of all points of the form (x, u) G A with x G {a, b}. 
The theory developed here is valid for more general boundary conditions. 

Let Pa,a,<5 denote the probability measure associated with the above con- 
tinuum percolation model on A. For a given configuration uj of deaths and 



Entanglement in the quantum Ising model 



19 



bridges on A, let k(uj) be the number of its clusters (subject to the top/bottom 
periodic boundary condition). Let q £ (0, oo), and define the 'continuum 
random-cluster' probability measure Pa,a,5,? by 

dV AjXM (u;) = ^q k ^d¥ AX8 (u;), lu £ n A , (5.2) 

for an appropriate 'partition function' Z. 

The theory of the continuum random-cluster model may be developed in 
very much the same way as that for the random-cluster model on a discrete 
lattice, see [12]. We shall assume the basic theory without labouring the 
calculations necessary for full rigorous proof. The details may be obtained by 
following minor variants of the relevant strategy for the discrete case. 

If fi is a probability measure and / a function on some measurable space, 
we denote by fi(f) the expectation of / under /i. 

The space Q A is a partially ordered space with order relation given by: 
Ui < u>2 if B(uji) C i?(cj 2 ) and D{ui) D D{uj2)- A random variable X : 
Q A — > R is said to be increasing if X(u) < X(u') whenever u < uj' . An event 
A £ T A is said to be increasing if its indicator function 1^ is increasing. Given 
two probability measures fii, fi2 on the measurable pair (Q^^a), we write 
Mi <st 1^2 if Mi (AT) < fi2{X) for all bounded increasing continuous random 
variables X : Q A — > R. 

The measures Pa,a,<5, 9 have certain properties of stochastic ordering as the 
parameters A, A, 5, q vary. There are two approaches to such stochastic 
inequalities, either by working on discrete graphs and passing to a spatial limit 
to obtain the continuum measures, or by working directly in the continuum. 
We shall not pursue this here, but refer the reader to [5] for a discussion of the 
case 5=1. The following two facts will be useful later. First, Pa,a,<5,9 satisfies 
a positive-association (FKG) inequality when q > 1, and secondly, 

W A ,x,s, q <st Pa,a,«5 when q > 1. (5.3) 

In the current paper we shall work mostly with finite-volume measures, 
that is, with measures defined on boxes of the form of A = [a, b] x [s,t]. We 
assume henceforth that q > 1. Having established the necessary estimates 
on such boxes, we will pass to the vertical infinite-volume limit as s — ► — oo, 
t — > oo. The existence of such a limit is not explored in detail here, but we 
note the following (see [T]). If we work on A with top/bottom wired or free 
boundary conditions, then the limit measures exist as a consequence of positive 
association (very much as in the lattice case, see [12]). Furthermore, the weak 
limit with top/bottom periodic boundary conditions exists and agrees with the 
first two limit measures whenever the latter are equal. A sufficient condition 
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for this is that the wired limit measure does not percolate. Since the limit of 
A as t — s^ooisa strip of bounded width, this condition is satisfied for all 
A, 5 £ (0,oo), and therefore the limit measures exist and do not depend on 
the choice of boundary condition. 

The situation is slightly less clear in the doubly- infinite- volume limit, as 
A f Z x I. The self-dual point for the continuum random-cluster measure 
on Z x K is given by X/5 = q, and thus one expects the free and wired limit 
measures to be equal at least whenever X/5 ^ q. It may be shown using 
duality that there is no percolation when X/5 < q, and it follows that the 
weak limits 

F XAq = lim ¥ AA5>q , q>l, 

exist if X/5 < q. We shall make no reference to this later. 

Just as the g-state Potts model may be coupled with a random-cluster 
model on a given graph, so may we consider a continuum Potts model on 
a box A = [a, b] x [s,t]. Let q £ {2,3,...}. We sample u according to 
^A.\,8,q, and we allocate a randomly chosen spin from the set {1,2, ... ,q} to 
each cluster of ou; the points of each cluster receive a given spin-state chosen 
uniformly at random from the q possible local states, and different clusters 
receive independent spin-states. We call the ensuing spin-configuration a q- 
state continuum Potts model, and a continuum Ising model when q — 2. When 
q = 2, by convention we take the local spin-space to be { — 1, +1} rather than 
{1,2}, and this is the case of interest in the current paper. We note in passing 
that the g-state continuum random-cluster model corresponds to a certain q- 
state quantum Potts model constructed in a manner similar to that of the 
quantum Ising model. 

The set of spin-configurations of the continuum g-state Potts model is the 
space Sa given as follows. Let T be the set of finite subsets of A. For D £ T ' , 
let J{D) be the set of maximal intervals of the time-lines that contain no point 
in D (subject to the top/bottom boundary condition on A). The space £a is 
defined as the union over D of the set of functions a : J(D) — > {1,2, ... ,q} 
with the property that <7( XiU _) 7^ &(x,u+) f° r a h (x,u) £ D. The corresponding 
probability measure on Sa is found by integrating over uj in the above recipe, 
as in the following summary. For a £ Sa, write D a for the set of points 
(x,u) £ A such that (J( X) u-) 7^ cr (x,u+)- The probability measure P associated 
with the continuum g-state Potts model on A is given by 

dW(a) = ^e XL ^dF 6 {D a ), a £ S A , 

where P5 is the law of an independent family of Poisson processes with intensity 
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5 on the time-lines indexed by [a, b], and 

is the total length of neighbouring time-lines where the spins are equal. Here, 
the summation is over all unordered pairs x, y of neighbours. We shall not 
develop the theory of such measures here, save for noting that P has the spatial 
Markov property (see [8j [H] for accounts of the spatial Markov property for 
a lattice model). For a G £a and a measurable subset S of A, we write as for 
the value of a restricted to S, and Qs for the cx-field generated by as- 

The above definition of the continuum random-cluster model is based on 
an assumption of free boundary conditions on left/right sides of the region 
A (we shall always assume top/bottom periodic conditions in this paper). 
More general boundary conditions may be introduced as follows. Let r be 
an admissible configuration of deaths and bridges off the box A. That is, r 
comprises a set D(t) of deaths and a set B(r) of bridges of (Z x R) \ A such 
that: the intersection of D(r) and B(r) with any bounded sub- interval of 
Z x R is finite, and no death is the endpoint of any bridge. For uj G Q\, we 
denote by (uj, t) the composite configuration comprising uj on A and r on its 
complement. We write ^\xsq f° r the continuum random-cluster measure on 
Q\ with the difference that the number k(uj) of clusters in (15.21) is replaced by 
the number k(uj, r) of clusters of (uj, r) that intersect A (subject, as usual, to 
the top/bottom periodic boundary condition). As in the lattice case, PX.A.a.g 
is stochastically increasing in r. One may consider also periodic boundary 
conditions. 

We extend this discussion now to boundary conditions defined in terms 
of spins rather than deaths/bridges. Let q > 2 be an integer. Let r be a 
boundary condition as above, and let n be a mapping from its clusters to the 
set {1,2, ... ,q}; that is, rj allocates a spin to each cluster of r, viewed as a 
configuration on (Z x R) \ A. Let the measure ^ v AX 5q ^ e gi ven as ^*A,\,s,q, 
conditioned on the event that no two points x,y G d h A with rj(x) ^ n(y) are 
connected. We now allocate spins to the clusters of the composite configura- 
tion (uj, t) by: if a cluster C contains a vertex y that is already labelled, the 
entire cluster of y takes that label, and if no such vertex exists, the spin of C 
is chosen uniformly at random from {1,2, ... ,q}, independently of the spins 
on other clusters. 
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6 Basic estimate for the slit box 

We consider next a variant of the above model in which the box A possesses 
a 'slit' at its centre. Let L > and Sl = [0, L] x {0}. We think of Sl as a 
collection of L+l vertices labelled in the obvious way as x — 0, 1, 2, . . . , L. For 
m > 2, P > 0, let A mj( a be the box [— m, m + L] x [— ^/?] subject to a 'slit' 
along Sl. That is, A m ( g is the usual box except in that each vertex x G Sl is 
replaced by two distinct vertices x + and x~. The vertex x + (respectively, x~) is 
attached to the half-line {x} x (0, oo) (respectively, the half-line {x} x (— oo, 0)); 
there is no direct connection between x + and x~ . Write Sf = {x ± : x G 
Sl} for the upper and lower sections of the slit Sl- We now construct the 
continuum random-cluster measure 4> m ,p on A mij g with top/bottom periodic 
boundary condition and parameters A, 5, q = 2. We shall abuse notation 
by using <$> m $ to denote also the coupling of the continuum random-cluster 
measure and the spin-configuration on A m / g obtained as above. An illustration 
of the slit box is presented in Figure [2J 

Let Q m ,/3 be the sample space of the continuum random-cluster model 
on A m p, and £ mjJ g the set of all possible spin-configurations. That is, £ m;/ a 
comprises all admissible allocations of spins to the clusters of configurations in 
VL m ^. For a G S mj/ g and x G Sl, write a^r for the spin-state of x ± . Let S L = 
{ — 1,+1} L+1 be the set of spin-configurations of the vectors {x + : x G Sl} 
and {x~ : x G Sl}, and write a\ = (<7+ : x G Sl) and cr^ = (cr~ : x G Sl). 

It may be checked from (14. ip that 

(pmA^L = e ~ 5 a l = e+ ) (e~|Pm(/3)|e + ), e",e + G S L , 

whence /(e + ,e~) = (j) m ,p{^L = e ~' °"l = e+ ) * s the function defined in (14. 3p . 
It is easily seen that a m ^, given in (14.51) . may be expressed as 

0>m,P = 4>mA a L = a l)- (6- 1 ) 

On recalling ([211]), by (Oil . 

MpLw - p L n (m) = ^^Sl _ (6 . 2) 

where c : — > C and 

^ = ^2 C ( a L) a L G 7^ L - 

The reduced ground state is obtained from p^iP) by taking the limit 
as (3 — > oo. By the remarks in Section [51 there exists a probability measure 
m such that 

4>m,f3 4>m as —> OO. 
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Furthermore, the are cylinder functions, and therefore, as (3 — > oo, 



m ,/3(c(a+)c(a L )) -> m (c(a+)c(a L )), (6.3) 

and 

«m,/3 ~» «m = 0m(o"i = (6.4) 



In order to prove Theorem 12.21 we seek the function c : — > C, with 

= 1, 




that maximises the modulus of (16. 2p . By splitting (16.21) into its real and 
imaginary parts, and applying the triangle inequality, we see that it suffices 
to consider functions c taking only non-negative real values. 

Here is the main estimate of this section, of which Theorem 12.21 is an 
immediate corollary with adapted values of the constants. 

Theorem 6.5. Let A, 5 G (0, oo) and write 9 = X/5. If 9 < 1, there exist 
a,C,M G (0,oo), depending on 9 only, such that the following holds. There 
exists 7 = 7(0) satisfying 7 > when 9 < 1 such that, for all L > 1 and 
M < m < n < 00, 



sup 



<f>m{c{(Tt)c(a L )) <p n {c{at)c{a L )) 



< CL a e~ im , (6.6) 



where the supremum is over all functions c : — > K. with \\c\\ = 1. The 
function 7 may be chosen to satisfy 7(6*) — >• 00 as 9 j 0. 

The condition < 1 is important in that it permits a comparison of the q = 
2 continuum random-cluster model on Z x R with the continuum percolation 
model. The claim of the theorem is presumably valid for 9 < 9 C where 9 C is 
the critical point of the former model. [It may be shown that 9 C > 2, and 
we conjecture that 9 C = 2, the self-dual point.] Similarly, Theorem 16.51 has a 
counterpart in d > 2 dimensions. 

We shall require for the purposes of comparison the following exponential- 
decay theorem for continuum percolation. Let A m denote the box [— m,m] 2 , 
and let / = {0} x [— |, |] be a unit 'time-segment' centred at the origin. 

Theorem 6.7. Let A, 5 G (0, 00). There exist C = C(X,5) G (0, 00) and 
7 = 7(A, 5) satisfying 7 > when A/5 < 1, such that: 

F x ,s(l <-> <9A m ) < Ce" 7m , m > 0. 

T/ie function j(X, 5) may be chosen to satisfy 7 — > 00 as 5 — ► 00 for fixed A. 
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Proof. Consider the continuum percolation process with parameters A, 5. The 
existence of such 7 is proved in [5] . That 7 — > 00 as 5 j (with A fixed) may be 
proved by bounding the cluster at the origin by a branching process. Consider 
an age-dependent branching process in which each particle lives for a length 
of time having the distribution of the sum of two independent exponentially- 
distributed random variables with parameter 5. During its lifetime, it has 
children in the manner of a Poisson process with parameter 2A, so that a 
typical family-size N has generating function 



The process is subcritical if E(N) < 1, which is to say that G' N (1) = AX/8 < 1. 
When this holds, the tail of the total number M of particles decays ex- 
ponentially, and similarly the aggregate lifetime U of the particles has an 
exponentially-decaying tail. See [131 E] f° r accounts of the theory of branch- 
ing processes. 

The branching process dominates C in the following sense. Identify the 
progenitor of the branching process and the origin of Z x 1. The length of 
the maximal death-free time-interval containing the origin has the distribution 
of the lifetime of 0. The number of bridges with an endpoint in this interval 
has the distribution of N. Each such bridge has endpoints of the form (0, s) 
and (x, s) where x = ±1. When we iterate this, we find that the number 
of bridges in the maximal death-free interval containing (x, s) is dominated 
(stochastically) by N. Arguing inductively, the number of bridges in the 
cluster C is dominated stochastically by the total size M of the branching 
process. 

The horizontal displacement of C is thus smaller (in distribution) than the 
total size M of the branching process. It is standard that the tail of M satisfies 
P(M > m) < Ce~ um for some C, v > depending on A, 5, and furthermore 
that v — > 00 if 5 I with A held fixed. The behaviour of v may be calculated 
exactly by elementary means, as follows. One may consider a variant of the 
branching process in which each particle has a lifetime with the exponential 
distribution, parameter 5, and has pairs of children at rate 2A while alive. The 
probability generating function of the total progeny may be found in closed 
form in the usual way (see [T3], Problem 5.12.11), and one obtains thus a 
sharp estimate for v via Markov's inequality. 

Similarly, the vertical displacement of C is smaller (in distribution) than 
the aggregate lifetime U of the particles in the branching process. Just as 
above, U has exponentially-decaying tail when E(N) < 1, and the constant 
in the exponent tends to infinity as 5 j for fixed A. 




s\ < 1. 
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Now, 

Pa,5(0 <-> dA m ) < P(M > m) + P([/ > m). 

A little more is needed for the theorem. The interval / is connected to a num- 
ber of bridges having the Poisson distribution with parameter 2A. The clusters 
generated by the ends of these bridges have sizes dominated (stochastically) 
as above, and the claim follows. □ 



In the proof of Theorem 16. 5[ we make use of the following two lemmas, 
which are proved in the next section using the method of 'ratio weak- mixing'. 

Lemma 6.8. Let A, 5 £ (0, oo) satisfy X/5 < 1. There exist constants 
a,Ci,C 2 £ (0, oo) such that: for all L > 0, m > 1, [3 > 2m + L, and all 
e + ,e" £ Sl, 

In the second lemma we allow a general boundary condition on A mj( g. 

Lemma 6.9. Let X, 8 £ (0, oo). There exist constants C, 7 £ (0, 00) satisfying 
< 7 < 1 when \/5 < 1 such that: for all L > 0, m > 1, (3 > 4(m + L + l), all 
events A C £ L x and all admissible random- cluster boundary- conditions 
t and spin boundary- conditions rj of A m ^, 



</>m,/?(OZ,ff L ) e A ) 



2 . 



< Ce 7^ m ; /or a = T,r], 



whenever the right side of the inequality is less than or equal to 1 . The function 
7 may be taken as that of Theorem \6. 7[ 

The above two lemmas are stated in terms of the box A m>/ g with top/bottom 
periodic boundary conditions. Their proofs are valid under other boundary 
conditions also, including free boundary conditions. We make use of this 
observation during the proofs that follow. 

The supremum in Theorem 16.51 may be handled by way of the next lemma. 

Lemma 6.10. Let \x be a probability measure on the finite set S . Let C be the 

class of functions c : S — > [0, 00) such that J2 s£S c(s) 2 = 1. Then 



ses V sgS 

with equality if and only if 



c( s ) = s £ S. 
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Proof of Lemma \6.1(A This is easily proved using a Lagrange multiplier. □ 



Proof of Theorem \6.5l Let < A < 5, and let 7 be as in Theorem 16.71 Let 
2 < m < n < 00 and take (3 > 4(m + L + 1). Later we shall let (3 — ■> 00. 
Since m>/ g < st n ,/3, we m &y couple m /3 and 4> n ^ via a probability measure v 
on pairs u^) of configurations on A n / j in such a way that z/(u;i < lo^ = 1. 
It is standard (as in [T2J [22]) that we may find v such that cl^ and lo 2 are 
identical configurations within the region of A m/ g that is not connected to 
<9 h A m „3 in the upper configuration ui 2 . Let D be the set of all pairs (cl>i,cl> 2 ) G 
f2 nj/ 3 x f2 nj/ g such that: uo 2 contains no path joining OB to <9 h A mjj g, where 
£? = [— r, r + L] x [— 2(r + L + 1), 2(r + L + 1)] and r (< \m) will be chosen 
later. We take free boundary conditions on B. The relevant regions are 
illustrated in Figure [3j 




Figure 3: The boxes A re ^, A m /3 , and S. 

Having constructed the measure v accordingly, we may now allocate spins 
to the clusters of lo\ and lo 2 in the manner described earlier. This may be done 
in such a way that, on the event D, the spin-configurations associated with 
uji and uj 2 within B are identical. We write 0\ (respectively, a 2 ) for the spin- 
configuration on the clusters of uj\ (respectively, cu 2 ), and af L for the spins of 
a-i on the slit Sl- 

For c : — > [0, 00) with ||c|| = 1, let 



so that 



.SI, 



c l^i 'lM°i,lJ 



*"m,f3 



I"n,l3 



l n,(3 



v(S c ;D)+v(S c ;D). 



(6.11) 



(6.12) 
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Here, D is the complement of D, and u(f; D) denotes u^fljj). 

We consider first the term v{S c ;D) in (16.121) . On the event D, we have 
that cr lL = <J^ L , so that 



HS C ;D)\ < 
By Lemmas 16.81 and I6.1CH 



a m,l3 



0>n,() 



Q"m,f3 



(6.13) 



< C 2 L> m ^(c(a+))0 mi/3 (c(^)) 



.14) 



where we have used reflection-symmetry in the horizontal axis at the interme- 
diate step. By Lemma [6.81 and reflection-symmetry again, 

a m ,/3 = E $™A a l = a L - e ) 



Therefore, 



0m,/3(c(>l)c(> L )) 
a m,/3 



(6.15) 



where C3 = C2/C1. 

We set A = {crj = o"^} in Lemma 16.91 to find that, for sufficiently large 
m > M'(A,<5), 



< Ce ~h m < -. 

By averaging over n, sampled according to (f) n> p, we deduce that 



>n,p\p L 



Ur = a, 



which is to say that 



a n,(3 . 



>"m,f3 



< Ce~^ m < -, 
2' 



(6.16) 
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We make a note for later use. By the remark after Lemma [6.91 inequality 
(16.151) holds also with 4>Tn,/3 replaced by the continuum random-cluster measure 
(f>B on the box B with free boundary conditions. Similarly, we may take C 
and M' above such that 



Wj 1 



a B 



< Ce~^ r <-, r> M'(X, 5), (6.17) 
2 



where a# = 0b(ctJ = er L ). 

Inequalities (I6.15P and (16 . 1 6[) may be combined as in (16.131) to obtain 

WS C ;D)\ < C 4 L 2a e~^ m (6.18) 

for an appropriate constant C 4 = C 4 (A, 8) and all m > M'(A, 5). 
We turn to the term v(S c \ D) in (16.121) . Evidently, 

\v(S c ; ~D)\ < A m + B n , (6.19) 

where 



0"m,/3 0>n,P 

There exist constants C5, M" depending on A, 5, such that, for m > r > M", 

a n,/3 

< ^Mc 5 B (c(a+ L ) C (a- L )) (6.20) 

CtB 

by Lemma [B31 with 4> m ,p replaced by 4>b, and (16.171) . At the middle step, we 
have used conditional expectation given the configuration r on A m> p \ B. By 
(16.151) applied to the measure 0b, there exists Cq = Cq(X, 5) such that 

— Mc(a+ L )c(a- L )) < C 6 L 2a . (6.21) 
a B 

Inequalities (I6.20p - (l6.2ip imply an upper bound for B n . 

A similar upper bound is valid for A m , on noting that the conditioning 
on D imparts certain information about the configuration 00% outside B but 
nothing further about U\ within B. Combining this with (I6.19I) - (I6.21I) . we 
find that, for r > M"'{\, 5) and some C 7 = C 7 (A, 5), 

HS C ;D)\ < u(D)C 7 L 2a . (6.22) 
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Let r = M'" to obtain by (15.31) and Theorem 16.71 that 

u(D) < C 8 {r + L)e~h m < C 9 Le~^ m , m > 2M'" , (6.23) 

for some C 8 , C 9 . We combine floUBJ) . (I6T22D . (I6T23|) as in (16321 . and let /3 -> oo 
to obtain ([63]) from (Q-([631), for m > max{M', M", 2M'"}. The constants 
C, 7 may be amended to obtain the required inequality. 

Finally, we remark that a, C, and M depend on A and 5. The left side 
of (16. 6p is invariant under re-scalings of the time-axes, that is, under the 
transformations (A, 5) i— > (Xrj, 5rj) for 77 G (0, 00). We may therefore work with 
the new values A' = 9, 5' = 1, with appropriate constants a(9, 1), C(6*, 1), 
M(0,1). □ 



7 Ratio weak- mixing 

Our proofs of Lemmas 16.81 and 16.91 make use of various couplings of random- 
cluster measures. Such couplings are fairly standard (see [121 122] for example) 
and have been utilised in [21 [3] in a study of ratio weak-mixing for random- 
cluster and spin models on discrete lattices. We follow in part the arguments 
of [21 E] , but we are not concerned here with the level of generality of those 
papers. 

Here is some notation. Let A be a box in Z x K (we shall later consider 
a box A with a slit Sl, for which the same definitions and results are valid). 
A path 7r of A is an alternating sequence of disjoint intervals (contained in 
A) and unit line-segments of the form [20,21], &12, [22,23], ^34, &2fc-i,2fc, 
[z2k, Z2k+i], where: each pair z 2 i, Z21+1 is on the same 'time-line' of A, and 
&2i-i,2i is a unit line-segment with endpoints z 2 i-i and z 2 i, perpendicular to 
the time-lines. Note that the equality z<n = Z21+1 is permitted. The path 7r 
is said to join zq and z 2 k+i- The length of ir is its one-dimensional Lebesgue 
measure, with ir viewed as a union of line-segments of M 2 ; note that each 
bridge of 7r contributes 1 to its length. A circuit D of A is a path except 
inasmuch as zq = z 2 fc+i- A set D is called linear if it is a disjoint union of 
paths and/or circuits. Let A, T be disjoint subsets of A. The linear set D 
is said to separate A and T if every path of A from A to T passes through 
D, and D is minimal with this property in that no strict subset of D has the 
property. 

Let to G f^A- An open path tt of uo is a path of A such that, in the notation 
above, the intervals [z2i,z 2 i+i] contain no death of u, and the line-segments 
&2i-i,2j are bridges of to. 

The (one-dimensional) Lebesgue measure of a measurable subset S of Z x IR. 
is denoted \S\. Let S and T be measurable subsets of A. The distance d(S, T) 
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from S to T is defined to be the infimum of the lengths of paths having 
one endpoint in S and one in T. Note that the distance function d depends 
on the choice of A (and, in particular, on the boundary conditions and the 
presence/absence of a slit). 

Let 0a denote the random-cluster measure on with parameters A, 5, 
q = 2 (with top/bottom periodic boundary condition). Let T be a measurable 
subset and A a finite subset of A such that A n T = 0. We shall prove a 
'ratio weak-mixing property' of the spin-configurations in A and T. In order to 
introduce the necessary couplings, we consider next a certain 'wired' boundary 
condition on A. Let denote the continuum random-cluster measure on A 
with parameters A, 5, q — 2, but subject to the difference that the set of 
clusters that intersect AUT count only 1 in all towards the cluster count k(u>) 
in (15. 2p . We call a 'wired random-cluster measure'. It is standard, just as in 
the discrete case, that may be used to generate a random spin-configuration 
on A corresponding to a continuum Ising model conditioned on having the 
same spin at all points in A LIT: let u be sampled according to 0, and allocate 
a randomly chosen spin from the spin set { — 1, +1} to each cluster of u, these 
spins being independent between clusters. 

Just as in the lattice case, one may use to obtain random-cluster measures 
with other boundary conditions. Let r G Er, and let Tj = {x G T : r(x) = i} 
for i = ±1. The corresponding random-cluster measure, denoted <p\ (as in 
Section [5]), is that obtained by: (i) the set of clusters intersecting V counts 
only 1 in all towards the cluster count in (15.21) . and (ii) we condition on the 
event that there exists no path joining T\ and T 2 . Since > st 0X, there exists 
a coupling k of the two measures with the property that k((ui,U2) '■ 0J\ > 
UJ2) = 1. It is natural to allocate spins to the clusters of u\ and u>2 in such a 
way that, whenever a cluster C of U2 is also a cluster of ui, and C D T = 0, 
then these two clusters have the same spin. 

One may carry out the above construction simultaneously for two (or more) 
r. Let t,t' G Er- We may find a coupling of 0, 0^, 0^ such that the first 
component is greater than each of the other two. That is, there exists a 
measure n on Q\ = {(cu, uj\, cu 2 )} such that: u> (respectively, u>i, u 2 ) has law 
(respectively, 0X, 0X), and k(lo > tui,^) = 1- 

Theorem 7.1 (Ratio weak-mixing). Let T C A be measurable, let A C A be 
finite such that An T — 0, and let D be a linear subset of A that separates A 
and T. Let A, 5 G (0, 00). For r, r' G Er and a G Ea, 




(7.2) 
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whenever the right side is less than or equal to 1, and where 

A <-►£), t 2 = ^(D~T). (7.3) 

The corresponding conclusion is valid when A is taken as the slit box A m 
Note in this case that the ti are given in terms of connection probabilities in 
the slit box. 

Proof. We adapt the methods of [2]. Let I (respectively, E) be the region of 
A reachable from A (respectively, V) along paths of A not intersecting D. 

Let r, t' G Sr and a G £a- We construct a coupling as follows, using the 
approach summarised prior to the statement of the theorem. Let uj have law 
(f). Let uj = uj t and uj' = uj t have laws cf) T k and <f)\ , respectively, and be such 
that uj,uj' < uj. Furthermore, we construct uj and uj' in such a way that, if 
uj G E2 = {D <r» T}, then uj, uj, and uj' are identical on D U /. 

To the clusters of uj, uj, uj' we assign spins in the usual manner, denoted 
a, a, a', respectively, such that: on the event E2, the functions a, a, a' are 
equal on D U /. For a reason that will be clearer later, we shall not work with 
the pair a , o' of configurations but instead with a pair p, p' defined as follows. 
First, we set 

Px = o~x, Px = °'x for x G -D U E. 

On the event F = {pn = pp}, we sample from the measure 0a given F to 
obtain a (random) configuration ( G S/, and we set 

Px = Px = Cx for x G /. 

On the complement of F, we sample p (respectively, p') according to the 
conditional law (f>\ given (p x : x G D U E) (respectively, (f>\ given (p' x : 
x G DUE)). By the spatial Markov property of the continuum Ising model 
alluded to after (15. 4ft . p (respectively, p') has law <fi\ (respectively, (p T A ), and 
furthermore: 

pj = p'j on the event {po = p'd}i (7-4) 

and 

n(p D = p' D ) = k(o- d = o-' D ) > k(E 2 ) = \-t\, (7.5) 

where k is the appropriate probability measure, and t 2 is as in (17.31) . 
Let H be an event satisfying 



H C {p A = p' A }. 



(7.6) 
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As in |2], if k(H) > 0, 

0a( cj a = a) k(pa = a) 



/<#n{p A = ct}) «(# | p A = a) 



| pa = a) k(H fl {p A = a}) 
k(H | p' A = a) 



(7.7) 



k(H I pa = a) 

It thus suffices, by an elementary argument, to prove that 

k(H I pa = «), k(H I p' A = a)<t (7.8) 

where 

t = *i + 2t 2 + * 1 + * 2 . (7.9) 

1 — Z\ — Zl2 

To see this, assume (17. 8ft with t<\. By (17.71) . 

i-t< ^A = q) < i 



1-t 



Now, 1/(1 - t) < 1 + 2t since t < |, and (O) follows. 

There are four steps in proving (17.81) . Let ^ (respectively, C^) be the 
cr-field generated by p^ (respectively, p' D ). Firstly, given that u E E\ = {A 
D}, the spin-vector op is (conditionally) independent of a A , whence 

\n(a D E A \ a A = a) — n(a D E A \ a A = a') | < t±, A E Q D , a' E E A , 

with t\ as in (17.31) . Averaging over a', we obtain 

k(<7d E A I cr A = a) — k(<j d E A)\ < ti, 

and hence, by the equidistribution of a and p, 

\k(p d E A I p A = a) - k( Pd E A)\ < t u AeQ d . (7.10) 

Secondly, let 

g = k(p d ^ p' D | Qd)i g = ^(pd ^ p' d I Q' D ), 

and, for a > 0, let H = H a be given as 

H a = {pD = p' D } n {g < a} n < a}, 
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where a will be chosen later. It is easily seen by (17.41) that H a satisfies (17.61) . 
By Markov's inequality and (17.51) . 

1 1 

*(g > a) < -n{g) < -tj, 
a a 

and therefore, since {g > a} G Qd, 

n(g > a \ pa = oi) < n(g > a) + 1\ by (I7.10p 

< + (7.11) 

a 



By a similar argument, 



K(g' > a I Pa = a) < ~t\ + t x . (7.12) 

a 



Thirdly, 



H 



^(pd ^ p'd, g < a I p A = «) < ess sup {n(p D ^ p' D | fo)l{ 9 <a}} 

= ess sup {gl{ g <a}} < a, (7.13) 

and similarly, 

K (pD ^ p'd, 9' < a I Pa = a) < a. (7.14) 

Finally, by (17^1) . 

{Pd = Pd> n {Pa = a} = {p D = p' D } n {p A = a}, (7.15) 
[this is where we use p, p' in place of a, a 1 }, and, by (IT. 12[) and (I7.14p - (l7.15p . 

«(po = p' D , g' > a I pa = a) < K>(g' > a \ p D = p Dl p' A = a) 

< 



n(g' > a 


Pa = a ) 


k(pd = p' D 


Pa = «) 



'■ + ^/° (7.16) 
1 — a — ti — 12/a 

On combining (17.111) . (I7.13p . (I7.16p . and setting a = t 2 , we obtain the first 
inequality of (17.81) with H = H a , and the second inequality holds similarly. □ 



Let A and T be disjoint finite subsets of A that are disjoint from d h A. 
Let D be an linear subset of A that separates A and V U d h A. Let a G Ea, 
(3, (3' G S r , and 77 G S 9 h A . By (O) applied to the sets A and V U <9 h A, 



^'Va = a) - rf'> A = a) I < 2<'> A = a), (7.17) 
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whenever t < ^ where 



t = h + 2t 2 + 



h + k 



1 — ti 



2U 



(7.18) 



and 



0(1} ^ru9 h A). 



(7.19) 



The suffix f3, r\ in (I7.17P indicates the composite boundary condition taking 
the values j3 on T and rj on <9 h A. We average (17.171) over (3' to obtain 



a A = a) - (j) A (a A — a) \ < 2t(j) v A (a A = a) 



(7.20) 



Now, 



^(cta = a) 



[a a = a, or = P) 
0>r = 0) ' 



Let A G Qa-i B G Qy be events with strictly positive probabilities. We 'multi- 
ply up' in (I7.20p and sum over a G A and (3 G -B to find that 



01(^)01(5) 



< 2t, 



Tj G S 



<9 h A; 



(7.21) 



whenever t < |. Upper bounds on t follow from the observation that <fi is 
stochastically dominated by the continuum percolation measure with param- 
eters A, 5 (cf. (15. 3p ). Equation (17.211) is a general statement of so-called ratio 
weak-mixing. 

By the same argument without the reference to the boundary <9 h A, 



<>k{a n B) 



4> A (A)(j) A (B] 



- 1 



< 2t 



A G Qa-, Be g r , 



(7.22) 



whenever t < |, where t is in d T . 1 8 [) with 



(7.23) 



and D is a linear set that separates A and T. 

The above ideas may be used to prove Lemmas 16.81 and I6.9[ for the first of 
which we argue as follows. Consider the box A m / g with slit Si,. Let K be an 
integer satisfying < K < ^L, and let A = {x + : x G Si, K < x < L — K} 
and T = {x~ : x G Si, K < x < L — K}. 
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Lemma 7.24. Let A, 5 G (0, oo). There exists C = C(X,5) G (0, oo) such 
that, if P > 2m + L, 



< Ce-^, e + e E A , e x G E r , 

whenever the right side is less than or equal to 1. The function j(X,S) may 
be taken as that in Theorem 6. 7 . 

The proofs are preceded by a type of 'finite-energy' inequality (see [2l [T2] ) . 

Lemma 7.25. Let S be a finite subset of A. For i 6 A \ S, f 6 Ej = 
{-1, and a G {-1, +1}, 

0aOs = e, o- x = a) > \<ph{a s = e)F A ^,s( x S )- (7.26) 

Proof. Let x G S, e G Eg, and a G { — 1,4-1}. Let E(e) be the decreasing 
event containing all uj G Qa such that: for all s, t G S, s t whenever e s ^ e t . 
Recalling the manner in which spins are associated with clusters, 

M°s = e) = A (2- fc(5) l s(e) ), e G E<j, (7.27) 

where k(S) is the number of clusters intersecting S. Similarly, 

M°s = e,a x = a)> A (2- fe(5+) l£ (e ) W), (7.28) 

where S + — S U {x}. Note that k(S + ) = k(S) + 1 when x ^ S. 
For any event A, 

M^ kiS+) U) = M^ HS+) )kA) = K$(A), (7.29) 

where K = <p A (2~ k ( s+ }) and is the continuum random-cluster measure on A 
with a wired boundary condition on S + , that is, all clusters intersecting S + 
are counted as one. By (17.291) and the FKG inequality applied to 0, 

M^ k{S+) ^E(e)l^s) = K${E(e) n{x^ S}) 

> K<j>(E(e))<j>(x S) 

= cp k {2-^ s+ h E{e) )i(x^s). 

Now k(S) < k(S + ) < k(S) + 1, so that, by (17271) -( 17251) . 

4>a(os = a x = a) > l4> A (as = e)$(x </» 5) 
and the claim follows by the stochastic inequality (15.31) . □ 
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Proof of Lemma 7. 24 Take D = {(x, 0) : x G [— m, 0) U (L, L + m}}, the union 



of the two horizontal line-segments that, when taken with the slit Sl, complete 
the 'equator' of A m ( g. Thus, D is a linear subset of A m>i g separating A and V. 
Since (p < st Pa,a,<5, by Theorem 16.71 there exist constants C, C depending on 
A and 5 alone, such that 

LL/2J 



and furthermore t% = t±. The claim now follows by f)7.22p . 



□ 



Proof of Lemma \6.8\ Let 7 be given as in Theorem 16.71 With K = [In L~\ , let 
o~f /K = (a^ . K < x < L — K). We may apply Lemma 17.251 as follows in order 
to compare the laws of the spin- vector and that of the reduced vector K . 
First, let x = (L,0), and let e + ,e~ G { — 1,+1} L+1 be possible spin-vectors of 
the sets S£ and S^, respectively. By Lemma [7.251 with S = S£ U S~[ \ {x + }, 



> 



0+ = e+ for y G S~£ \ {x + }, a' L = e -)P Am „ ^ S). 



Now, Pa 77I i3,\,&{ x w *S*) is a t least as large as the probability that the first event 
(death or bridge) encountered on moving northwards from x is a death. That 
is, 

On iterating the above argument, we obtain that 

<PmA a L = e+ > °L = O > ( 0/0 x I n ) ^mA a iK = e ^ 



2(2A + S) 



Ki a L,K 



(7.30) 

where e^- is the vector obtained from by removing the entries labelled by 
vertices x satisfying < x < K and L — K < x < L. In summary, there exist 
C, a G (0, 00) depending on A, 5 such that, for e ± G Si, 

CL~ 2a c 



A a L,K - 6 Ki a L,K 



e K ) < (pmA <J L 



cr, 



o 

< <f>mA a tK = e K' a LJ< = e 



Set A = {x + : x e S L , K < x < L - K}, T = {x~ : x G S L , K < x < 
L — K}, and apply Lemma [7.241 to obtain that there exists C = C(X, 5) < 00 
such that 



<PmA a 



L,K 



€ Ki a L,K 



-K) 



iA a L,K - e K) ( PniA a L,K 



~K> 
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whenever (say) the right side is less than or equal to |, say for L > L (A, S). 
By Lemma [7.251 again, for suitable C, a, 



C'lT 



+ 



± 



P\ a L,K - e K 



< 



e ± )< 



-K) 



The claim now follows for L > L , with suitable values of C±, C2, a. We may 
adjust the constants to obtain the required inequality for all L > 0. □ 

Proof of Lemma \6. 9[ Let A = S~£ U S£ and T = d h A m ^. Let k = |m and 
assume for simplicity that k is an integer. [If either m is small or k is non- 
integral, the constant C may be adjusted accordingly] Let D be the circuit 
illustrated in Figure HI comprising a path in the upper half-plane from (— k, 0) 
to (L + k, 0) together with its reflection in the x-axis. 



By Theorem 17.11 



(eV~)) 



(f>mA( a L> a L. 



O) 



< 2t, 



a — rj,T, e G Ej 



whenever t < |, with t as in (17.181) . We 'multiply up' and sum over (e + , e ) £ 
A to obtain 

^(o-a e A) 



>A e A) 



< 2t, 



(7.31) 



whenever t < |. 



By < st Pa,a,5- Let /? > 4(m + L + 1). It is a straightforward 

consequence of Theorem 16.71 that there exist C,C',c' > 0, depending on A, 5 
only, such that 



LL/2J [L/2] 
i=0 i=0 

and similarly, 

\k+L/2\ 
i=0 

with 7 given as in Theorem 16.71 The claim of the lemma follows. 



(7.32) 



(7.33) 
□ 



8 Disordered interactions 



We have so far assumed that the spin-couplings \ x , x +\ and the field-strengths 
8 X appearing in the Hamiltonian (11.21) are constant. The situation is more 
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2m + L 



Figure 4: The circuit D is approximately a parallelogram with A at its centre. 
The sides comprise vertical steps of height 2 followed by horizontal steps of 
length 1. The horizontal and vertical diagonals have lengths 2k+L and (order) 
4k + 2L respectively, where k = |m. 



complicated if: either the environment of couplings and strengths vary about 
the space Z, or they are random (in which case the model is said to be dis- 
ordered). The arguments of this paper may be applied in each case, and the 
outcomes are summarised in this section. 

Suppose first that the A^+i and 5 X are non-constant. The fundamental 
bound of Theorem 16.51 depends only on the ratio 9 — X/5, and the connection 
probabilities of the continuum random-cluster model are increasing in the 
X x ,x+i and decreasing in the 5 X . One may therefore check that the conclusions 
of the paper are valid with 7 = 7(A, 5) whenever 

K, y /S x < X/5, y = x - 1, x + 1, x G Z. (8.1) 



Hence, in the disordered case where ( 18. ip holds with probability one, the 
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corresponding conclusion is valid. 

We turn to the situation in which (18.11) does not hold with probability 
one. Suppose that the A xx +i, x G Z, are independent, identically distributed 
random variables, and similarly the 5 Z , z G Z, and assume that the X x<y 
are independent of the 5 Z . We write P for the corresponding probability 
measure, viewed as the measure governing the 'random environment', and A, 
A for a typical spin-correlation and field-strength, respectively. Conditional 
on A = (X x , x +i : x G Z) and 5 = (5 Z : z G Z), we write Pa,5 for the probability 
measure of the associated continuum percolation process. In applying the 
methods of this paper within the random environment, one needs to deal with 
sub-domains of Z where the environment is not propitious for the bound of 
Theorem 16.51 As before, we perform a comparison of the continuum random- 
cluster model and continuum percolation in a random environment, and we 
shall appeal to the following theorem of [T7] (see also Theorem 1.6 of pQ). 

For (x, s), (y, t) G Z x R and q > 1, let 

d q (x, s;y,t) = max{ \x — y\, (ln + \s — t\) q }, 
where ln + x = maxjlnx, 0}. 

Theorem 8.2. | I17| Consider continuum percolation on Z x R in a random 
environment satisfying 

T = max {P([ln(l + A)f) , P([ln(l + A" 1 )]' 3 ) } < oo, 

for some 

(3>5 + \\f2. (8.3) 

There exists Q = Q{(3) > 1 such that the following holds. For q G [1, Q) and 
7 > 0, there exists e = e(/3, T, 7, q) > and r\ = r](/3, q) > 1 such that: if 

P([ln(l + (A/A))f) <e, (8.4) 

there exist identically distributed, positive random variables D x G L r? (P), x G 
Z, such that 

Pa,5 {(x, s) <r-> (y, t)) < exp [--fd q (x, s; y, t)] if d q (x, s; y, t) > D x , (8.5) 
for (x, s), (y,t) G Z x R. 

The lower bound (18. 3p for (3 is enough to imply that P(D^) < 00 for some 
rj > 1. The larger /?, the larger 77 may be taken. 

For the remainder of this section we assume that the conditions of the 
above theorem are valid, and we shall work with the conclusion ( 18.51) . with 
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q — 1 , 7 > 1 , and the D x given accordingly. We let L > 8 and K = [In L] . 
and consider the event 



L-K 



A L = P| {Ac < min{x, L - x}}, 



i=.fsr 

noting that 



P(A £ ) > 1 - 2 ^ P(D > 



_ x) 

x=K 

where D has the distribution of the D x . Since P{D) < oo, 

P(A L ) -»■ 1 as L -> oo. (8.6) 

An estimate for the rate of convergence may be obtained (here and later) by 
the fact that P(D V ) < oo for some rj > 1. 

We comment next on the adaptation of our earlier results to the disordered 
setting. Theorem 17. II holds within the random environment, without change. 
The conclusion of Lemma 17.241 is valid with K = [In L] whenever the event 
Al occurs. Lemma 17.251 holds unconditionally. The conclusion of Lemma 16.81 
holds on Al with the lower bound C\L~ a replaced by CXl and the upper 
bound C2L 01 replaced by (CXl) -1 , with C a constant and 



q S x + X x ,x-1 + ^x,x+l 

where, in the notation of the proof of Lemma 16.81 6 = (S£ \ A) U (S^ \ T). 
Now, 

K-l L 
x=0 x=L-K+l 

where 

. ^x,x—l ^x ,x+l 

In H 



8 X 

The two summations in (18. 7p are independent of one another, and each is the 
sum of a 1-dependent sequence of random variables. Also, 

Z x <ln (l + ^pl) + In (l + ^f^ 
so that, by (18.41) and the Minkowski inequality, 



y/PW) < 2 JP([M1 + (A/A))] 2 ) < 00. 
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By the central limit theorem for 1-dependent sequences (see, for example, 
Theorem 19.2.1 of p]), 

P{B P L ) -»• 1 as L -> oo, (8.8) 
where = {Xl > L~ p } and p G (0, oo) satisfies 

p > AP(Z ). (8.9) 



Some changes are necessary to the proof of Lemma 16.91 reflecting the fact 
that the decay in (18. 5p is sub-exponential in time. The circuit illustrated 
in Figure H] is generated by translation, discretisation, and reflection of the 
Cartesian line y = 2x. In the disordered setting, we work instead with the 
curve y = e x , and we assume j3 > 5e m+ 2 L . We define two further events 
that depend on the environment. Assume for simplicity that m is even, write 
k = hm, and let 



2""' 



CL, m = f]{ D x < |min{A; + x, L + k - x}}, 

2=0 

L+k 

P > L,m — f~l {Dx < min{m + x, L + m — x}}. 



x=—k 

In the current setting, (I7.32p becomes 

ti < Cie _ 4T m on the event Cl,™,, 
for some constant C\ depending on 7. Similarly, (I7.33P is replaced by 

t\ < C 2 e~^ m on the event D L ^ m . 

An amended version of Lemma l6\9l thus holds, so long as the event CL, m ^P ) L,m 
occurs. 

We estimate P{CL, m H P>l,w) as follows. First, since P(D) < 00, 

P(C L , m ) > 1 -2^P{D X > l(k + x)) -> 1 as m ^ 00. (8.10) 

x=0 

Similarly, 

P(D L , m ) > 1 - 2 P{D X >m + x)^l as m -> 00. (8.11) 
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Suppose that Al fl B p l fl Cl,™ H -Dl,™ occurs for some p satisfying ( 18. 9ft . 
The principal estimate (12.31) follows with CL a replaced by CL P as above. On 
the above event, the proof of Theorem 12.81 may be followed to obtain the 
logarithmic decay of entanglement. Note from (18.61) and (18. 8p that P(Al fl 
B P L ) -»• 1 as L -> oo, and by (IHTTOD - flHTTTD that P(C L>m n L> L)TO ) -> 1 as 
m — > oo. 

Proof of Theorem \8.2[ This is essentially Theorem 1.1 of [17] with d = 1, 
subject to two differences: the right side of (I8.5P is expressed differently in 
[T7] . and the condition on (3 is different. The present statement is obtained 
as follows from the proof of (T7J, using the notation of that proof. With (3 
satisfying (18.31) and a = 1 + y/2, we pick p > 2a and v = q^ 1 satisfying (3.3) 
of [TTJ . Let K x denote the minimal k\ in the second paragraph of the proof of 
Theorem 3.3 of [17]. As there, 

P(K X > r) < r > 1, 

where c is a constant, and (L r : r > 1) is a sequence of positive reals given 
by L r = L aT for some large L. Let D x = oLk x - Inequality (18.51) holds by the 
argument of [T7J. Furthermore, for rj > 1, P(D^) has the same order as 



oo oo 

W ^P(L Kx >a)< W L^' +1 ■ jj—^, (8.12) 

a=0 r=0 

which is finite whenever 77 — 1 is small and positive. □ 
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